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Electronic structure of bcc Np, fee Pu, Am, and Cm pure metals under pressure has been investigated 
employing the LDA+t/ method with spin-orbit coupling (LDA+f/+SO). Magnetic state of the actinide ions 
was analyzed in both LS and jj coupling schemes to reveal the applicability of corresponding coupling bases. 
It was demonstrated that whereas Pu and Am are well described within the jj coupling scheme, Np and Cm 
can be described appropriately neither in {ma}, nor in {jrrij} basis, due to intermediate coupling scheme 
realizing in these metals that requires some finer treatment. The LDA+[/+SO results for the considered 
transuranium metals reveal bands broadening and gradual 5/ electron delocalization under pressure. 



PACS: 71.20.-b, 71.27.+a, 71.50.-y 

Prominent structural transition from S- to a-phase 
,of plutonium gives substantial total volume contrac- 
tion [T]. In americium and curium the similar transi- 
tions were also found under high pressure or increasing 
temperature |2l [3j . These volume collapses are usually 
related to the drastic electron delocalization. In curium, 
additional magnetic stabilization of intermediate phases 
■was also found [31 H] . 

In many strongly correlated materials Coulomb and 
Hund interactions are dominant, whereas spin-orbit 
(SO) coupling is comparably smaller and usually can be 
'treated within various perturbative techniques. In the 
lease of 5/ electronic shell of actinide metals all these 
'three terms in Hamiltonian are of comparable strength. 
That results in a very sensitive (non)magnetic ground 
state, but also in coupling schemes varying from usual 
LS (Russel-Saunders) one. An intermediate or jj cou- 
pling schemes were found to be more appropriate for 5/ 
electrons in transuranium metals O [6l [7l [8] . 

In actinide elements spin-orbit coupling is stronger 
than exchange Hund interaction, and hence the jj cou- 
pling scheme can be valid with a well defined total 
moment J, but in this case spin S and orbital L mo- 
ments are not well defined. Then the basis of eigen- 
functions of total moment operator {jrrij} is the best 
choice, since the matrix of spin-orbit coupling operator 
and occupation matrix are diagonal in this basis. How- 
ever the exchange interaction (spin-polarization) term 
in the Hamiltonian is not diagonal. 



In some 5/ elements intermediate coupling scheme 
is realized, then the occupation matrix has nondiagonal 
elements in both {ma} and {jmj} orbital bases. There- 
fore, both terms in the Hamiltonian: spin-orbit coupling 
and exchange interaction, should be taken in a general 
nondiagonal matrix form. 

Numerous band methods and approximations have 
been applied to describe magnetic and spectral proper- 
ties of transuranium metals, see for review. Nonmag- 
netic ground state of pure plutonium metal observed ex- 
perimentally [S] was reproduced in the electronic struc- 
ture calculations in the LDA-f C/-I-SO calculations [TU] 
(local density approximation supplemented with the 
Hubbard f/-correction and spin-orbit coupling). In 
these calculations the exchange interaction was found 
to be the reason of artificial antiferromagnetic ordering 
in various LSDA+C/ investigations. Also non-magnetic 
ground state was found in around-mean-field version of 
the LDA+U method [11], and later on within LDA -|- 
Hubbard I approximation [T2], and also hybrid density 
functionals with a dominant contribution of HF func- 
tional [T3'. Recently, the reliability of these results 
was supported by more detailed analysis of exchange 
interaction 14J, and also in the LDA-fDMFT calcu- 
lations [15 , supplementing LDA with the Dynamical 
Mean-Field Theory [16] (DMFT) [ITl O O [20] . 

While consistent interpretation of spectroscopic data 
is not found yet |21j, the LDA+C/+SO method provides 
realistic magnetic state and allows one to estimate elec- 
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FIG. 1. The density of states: total (solid line) and 5/ 
(shaded area) states of bcc Np, fee Pu, Am, and Cm 
from the LDA+f/+SO calculation. The Fermi energy 
corresponds to zero. 



trical resistivity in actinide metals [22] and alloys [25] 
in good agreement with experimental data. 

In this paper electronic structure of bcc Np, fee Pu, 
Am, and Cm pure metals was calculated within the 
LDA+C/+SO method introduced in detail in Ref. [TO] . 
In the method the exchange interaction (spin polariza- 
tion) term in the Hamiltonian is implemented in a gen- 
eral nondiagonal matrix form regarding the spin vari- 
ables. This form is necessary for correct description of 
5/ electrons behavior for the cases of jj and intermedi- 
ate coupling types. 

In the LDA+J7+S0 method accounting for strong 
electron correlations includes Coulomb interaction ma- 
trix elements which can be expressed via direct U and 
exchange Jh Coulomb parameters. A reliable way to 
estimate these parameters is provided by the constrain 
LDA approach [25] . To calculate a value of U parame- 
ter, in this procedure one evaluates a screened Coulomb 
interaction of 5/ electrons which requires the choice of 



screening channels taken into account. For s, p, and d 
channels for bcc Np and fee Pu, Am, and Cm the con- 
strain LDA calculations resulted in the Coulomb param- 
eter value U — A eV [THl [13 • The exchange Coulomb 
parameter Jh is evaluated as the difference of inter- 
action energy for the electrons pairs with the opposite 
and the same spin directions. Parameter Jh does not 
depend on the screening channels choice. For neptu- 
nium and plutonium, the value of Hund exchange pa- 
rameter Jh was calculated as Jh = 0.48 eV [TU], for 
americium as 0.49 eV, and for curium as 0.52 eV [22]. 
In this work in all LDA+C/+SO calculations we used 
these exact values of Jh for each actinide metal, since 
(non)magnetic state of actinide metals is sensitive to the 
value of Jh [ID] . 

Strong spin-orbit coupling of 5/ electrons results in 
a splitting of the / band into subbands corresponding 
to the values of total moment j=5/2 and j=7/2. The 
value of this splitting is 1 - 1.5 eV. Taking into account 
Coulomb correlations via the LDA+U correction does 
not change qualitatively the band structure, only the 
separation between subbands increases from 1.5 eV to 
5 eV according to the value of C/ = 4 eV. 

The density of states (DOS) from the LDA+[/+SO 
calculations for Np, Pu, Am, and Cm at ambient pres- 
sure are shown in Fig. [Tj In all DOSs one can distin- 
guish two groups of bands: j—5/2 at the lower energies 
and j=7/2 at higher ones. In Np the j—5/2 subband 
is partially filled. From Pu to Cm the Fermi level is 
shifted upward from the upper slope of j=5/2 subband 
in Pu and crosses the j=7/2 subband in Cm due to 
the increasing number of / electrons. These results are 
for the metals at ambient pressure in the cubic phases. 
We used the same lattice parameters as in our previous 
work [25] . 

To model the pressure, we assume it to be uniform 
and applied as a cell volume contraction. While exper- 
imentally crystal structure of these metals transforms 
from cubic to complicated structures like orthorhombic 
and monoclinic, for the electronic structure and mag- 
netic state of the actinide ion the main result of the 
applied pressure comes from the contraction of the unit 
cell volume per ion. In Figs. |2]-|5] one can find density of 
states (DOS) for the cubic phases under pressure, real 
phases with corresponding volumes per ion are referred 
to in brackets. In the actinide metals under investiga- 
tion the cubic phases do not correspond to the largest 
volume per ion, for this reason we also consider some 
volumes exceeding the cubic volume at ambient pres- 
sure. 

Neptunium metal. - Magnetic ground state for Np 
at ambient pressure was found in our calculation (see 
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FIG. 2. The partial DOS for the j=5/2 (shaded area) 
and j—7/2 (dashed hne) subbands of bee neptunium 
under pressure obtained from the LDA+f/+SO calcu- 
lations. For each curve the relative volume V/Vo is 
given. The Fermi energy corresponds to zero. 



Table 1) . Large values of the off-diagonal elements (OD) 
in Table 1 in both {mcr} and {jm.j} bases evidence for 
intermediate coupling type in this actinide metal. Four 
largest eigenvalues close to the unit give /^ configura- 
tion. Upon applied pressure the density of states at the 
Fermi energy increases, see Fig. ^ effective magnetic 
moment fJ.'^'ft ranges from 2.55 fiB to 2.33 fiB- 

Plutonium metal. - The LDA+C/+SO calculations 
for metallic Pu fee structure (fee phase in Pu is named 
(5-phase) at all volumes gave a nonmagnetic ground state 
with zero values of spin S, orbital L, and total J mo- 
ments [lOj in agreement with numerous experimental 
data [S]. The occupation matrix has six eigenvalues 
close to unit, see Table 1 (for details, see Ref. [l^), and 
is nearly diagonal in the {jruj } basis of eigenfunctions 
of total moment J. That gives a /^ configuration of Pu 
ion in jj coupling scheme. In Fig.[3]the partial densities 
of states for /^/^ and /^^^ subshell of Pu in all volumes 
(applied pressures) are presented. The DOSes as well as 
the occupation matrix demonstrate almost completely 
filled /^/^ band with the Fermi level on the top of it 
and an empty /^" band. The separation between the 
centers of these bands is « 5.2 - 5.4 eV. For small vol- 
umes j=5/2 subband does not change its position but 
its bandwidth becomes approx. 1 eV larger. 
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FIG. 3. The partial DOS for the j = 5/2 (shaded area) 
and j = 7/2 (dashed line) subbands of fee plutonium 
under pressure obtained from the LDA-|-(7+S0 calcu- 
lations. For each curve the relative volume V/Vo is 
given. Real phases with corresponding volumes per ion 
are cited in brackets. The Fermi energy corresponds to 
zero. 

Americium metal. ~ The calculates DOSes of Am 
are shown in Fig. H] One can see that j=5/2 subband 
is fully occupied that corresponds to /^ configuration 
whereas j=7/2 subband is almost empty. The occu- 
pation matrix has six largest eigenvalues of close to 
unit, see Table 1. The Fermi level in Am is shifted 
towards j=7/2 subband comparing to Pu due to addi- 
tional valence electron (see Fig. H]) that occupies s, p, 
and d states (not shown separately). Note, that the 
5/ bandwidth in Am is smaller than in Pu due to the 
larger cell volume per ion. Having a delicate balance 
between SO and exchange interactions, calculated 5/- 
DOS occupied 5/ band is centered around 4 eV in agree- 
ment with Am photoemission spectra that demonstrate 
large density of states in the range (-4 eV; -2 eV) [^S] . 
The LDA+C/+SO calculations for Am gave a nonmag- 
netic ground state in all volumes with the 5/ shell with 
S = L = J = Om agreement with experimental data [5] . 

Curium m,etal. - In curium /^ " subband is partially 
filled, providing a peculiar empty state well below the 
Fermi level, see Fig. \5\ The LDA-f C/-I-SO calculation 
results in S and L moments listed in Table 1. The value 
of effective magnetic moment fJ.'^'ff calculated from the 
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TABLE 1. Electronic configuration of 5/ shell in Np, Pu, Am, and Cm ions in the cubic phases calculated within 
the LDA+[/+SO method. The largest values of occupation matrices off-diagonal elements ODls and ODjmj in the 
corresponding basis sets are given in the second and third columns. Then the seven largest eigenvalues of occupation 
matrix are presented. The columns from the eleventh to thirteenth show the calculated values for spin (5*), orbital (L), 
total (J) moments. 



Metal OD 



LS 



OD, 



Largest eigenvalues 



S 



L 



Np 


0.36 


0.46 


Pu 


0.45 


0.01 


Am 


0.47 


0.02 


Cm 


0.31 


0.45 



0.05 0.09 0.26 0.89 0.91 0.91 0.92 
0.03 0.92 0.92 0.93 0.93 0.93 0.93 
0.07 0.97 0.98 0.98 0.98 0.99 0.99 
0.10 0.99 1.00 1.00 1.00 1.00 1.00 



J 



1.40 


4.68 


3.28 
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FIG. 4. The partial DOS for the j = 5/2 (shaded area) 
and j — 7/2 (dashed line) subbands of fee americium 
obtained from the LDA-|-[/-|-SO calculations. For each 
curve the relative volume V/Vo is given. Real phases 
with corresponding volumes per ion are cited in brack- 
ets. The Fermi energy corresponds to zero. 
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FIG. 5. The partial DOS for the j = 5/2 (shaded area) 
and j — 7/2 (dashed line) subbands of fee curium ob- 
tained from the LDA+[/+SO calculations. For each 
curve the relative volume V/Vo, where Vo is the vol- 
ume at ambient pressure, is given. Real phases with 
corresponding volumes per ion are cited in brackets. 
The Fermi energy corresponds to zero. 



total moment value J by the method used in Ref. [10] 
for different volumes gradually decreases from 7.44 /ib 
for the largest volume (corresponding to the volume per 
ion in the Cml phase) to 6.86 fiB (CmlV). This result 
corresponds to the model calculations assuming in pure 
LS coupling the magnetic moment of the curium ion to 
be 7.94 ^B, but for a reahstic model of intermediate cou- 
pling lowering up to 7.6 hb [1S]j while experimentally 
the magnetic moment was reported as 7.85 /xb [55]. 

In conclusion, we have calculated the electronic 
structure of bcc Np, fee Pu, Am, and Cm pure metals 
within the LDA-f C/+SO method under applied pressure 
in assumption of uniform unit cell contraction. Static 



mean-field approach used in this work results in gradual 
(featureless) electron delocalization and bands broaden- 
ing under pressure. While jj coupling scheme is well 
suited for description of Pu and Am ground state, for 
Np metal only intermediate coupling scheme seems to 
be appropriate, whereas in Cm the intermediate cou- 
pling is closer to LS-type. We suggest accounting for 
dynamical correlations effects will allow one to repro- 
duce not only the delocalization of 5/ electrons under 
applied pressure but also will resolve some fine feature 
of this process. 
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